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Abstract. The well-known analytical solution of Burgers' equation is extended to curvilinear coordinate 
systems in three-dimensions by a method which is much simpler and more suitable to practical applications 
than that previously used [22]. The results obtained are applied to incompressible flow with cylindrical 
symmetry, and also to the decay of an initially linearly increasing wind. 
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1. Introduction 

On 

Burgers' equation is a well-known example of a non-linear partial differential equation whose solution can 
be constructed from a linear partial differential equation. It is, to the best of our knowledge, the only such 
example. Much of the interest in this equation arises because it is a very simple form of the Navier-Stokes 
equation in the one-dimensional, cartesian, time-dependent, compressible, viscous limit [17]. An early study 
derived two steady-state solutions [1], and then the equation was used in studying the decay of free turbulence 
[4,5]. The importance of the equation is due to the non-linear term, wfj, which allows the calculation of 

CO ' 

the modification of the velocity due to the exchange of momentum between a variety of different length 
scales. Burger's equation is the simplest type that shows the complicated interplay between the non-linear 
steepening and diffusion of a wave. 

The first full solution of the one-dimensional Burgers' equation was found independently by both Cole [7] 
and Hopf [12]. Recently, the solution has been extended to n-dimensions by using group action on coset 
bundles [22]. The basic idea of this approach is to determine coset spaces for a chosen group and for a 
chain of closed subgroups, and then to construct a bundle out of these spaces. Next, the action of the group 
on the bundle is prescribed and a class of tensor-valued functions called invariants on the coset bundle is 
introduced. Because these group invariants take the form of differential equations, the method has been 
used to obtain a tensor version of the n-dimensional Burgers' equation and to show that its solutions can be 
constructed from a linear tensor diffusion equation. It has also been demonstrated that the procedure gives 



a nonlinear tensor equation for a generalized Burgers' equation with a higher than second order derivative. 
In this paper, we examine extensions of the basic class of solutions of the vector Burgers' equation, which we 
discovered before we became aware of the work of Wolf et al., to more than one dimension and to curvilinear 
coordinate systems. The main advantage of our approach is that it is extremely simple from a mathematical 
point of view and that the results we obtain have the potential to be used in a number of applications (see 
Section III for two simple applications). It must be added, however, that the solutions obtained in this paper 
represent a subclass of more general solutions to the tensor Burgers' equation discussed by Wolf et al. [22]. 

In what follows, we recapitulate part of the historical derivation presented first by Cole [7] and Hopf [12]. 
This derivation is central to our results in Sections II and III. 

The equation to be solved is 

du du d 2 u 

where v is a physical constant (viscosity). A successful transformation uses 



u=|^, (f> = <l>(x,t). (2) 



Substituting into eq. (1) and integrating term- by-term with respect to x 

d(j) 1 fdc/)\ 2 d 2 4> 
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where the function of integration is omitted. Because eq. (3) is invariant under the transformation 

x ^ ax, (4a) 

t -> a 2 i, (46) 
(a=constant), as is the diffusion equation, this suggests that there exists a solution of the form 

</>(x,t) = <j>(0(x,t)), (5) 
where 6 is a solution of the diffusion equation 

de d 2 e 

dt = "a?" (6) 

Substitution into eq. (3) while subtracting eq. (6) gives 

UJ ^ (7o) 

whose solution is 

4>{6) = -2v\n{6 - a) + c 2 , (76) 



so that by choosing ci=0 

<^) = —-. (8) 

Quite an interesting literature has developed around applications of the solutions in pure mathematics [20] , 
magnetohydrodynamics [9], astrophysics [14,16], and cosmology [10,15]. Lighthill [18] extended the basic 
transformation to the case where the x coordinate is measured in the frame of reference moving in the same 
direction as the wave at the undisturbed speed of sound. A valuable book on Burgers' equation and related 
topics was recently published but is, unfortunately, already out of print [11]. 

In the following, wc formally present the derivation of the solutions to the vector Burgers' equation by using 
a generalized Colc-Hopf transformation. 

2. The Vector Extension of Burgers' Equation 

We begin with the vector equivalent of Burgers' equation 

— + u • Vw = vV 2 u. (9) 
at 

In order to evaluate the vector Laplacian we must use 

V 2 u = V(V • u) - V x V x u, (10a) 

and for the inertial term 

u 2 

u-Vu = Vy - u x V x u. (106) 
Then eq. (9) greatly simplifies as long as 

uxVxu = iA7xVxu. (10c) 

This appears to be an equation that represents quite general flows, but the only solutions are irrotational. 
This can be seen by taking the divergence of both sides of eq. (10c) and using 

V ■ (Ax B) = B -V x A- A-X7 x B. (lOd) 

In what follows we will assume that the flow field is irrotational as this is necessary to simplify eq. (9) 
and to derive the extension of the Cole-Hopf transformation. The assumption is satisfied trivially for the 
one-dimensional Cartesian problem. 

Vxu = 0, (11a) 

so that 

u = V(p. (lib) 

Eq. (lib) is the extension of eq. (2) to three dimensions and, together with eqs. (10) allows the simplification 
of eq. (9) to 



dej) (V0) 2 „ 2 



0. (12) 



The meaning of eq. (12) is that the quantity enclosed by brackets must be a function of time, only because 
its gradient is zero in all space so that 

The form of the unknown function E(t) is irrelevant for our discussion as it cannot affect the velocity. We 
simply define a new potential 

01 = <f> - J E(t)dt, 
where the gradients of <f) and 4>\ are identical. 

Before introducing the Cole-Hopf assumption, it is useful to state the following vector identity: 

V« - §V0, (14a) 

which is true for any orthogonal, curvilinear coordinate system. Taking the divergence of both sides with 
a = (f> and using cq. (14a) again with a = ^ yields 

^ = W {v6)2 + fe v2e - (146) 



We will also need the diffusion equation 



(15) 



Now we apply eqs. (5), (14b), and (15) to eq. (13) to generate eq. (7a), again, so that 

u = - T ve. (16) 

This clearly shows that the Cole-Hopf transformation goes through for quite general coordinate systems. 
Cole [7] originally discussed the vector solution for the case of cartesian coordinates and he indicated that it 
applied to some special solutions in higher dimensions. 

The usefullncss of eq. (16) can be shown by exploiting the simple symmetry of the following example. 

3. Cylindrical Symmetry 

Using only a radial velocity, Burgers' equation in cylindrical coordinates with axisymmetry may be written 
as 

fin fin \ 1 f) I fin \ it 

(17) 



du du 

Equation (17) can be written as 



dt ^ U dr 



r dr \ T dr 



du du d ( 1 d \ 

Tt +U Wr =y d-r\rd-r {ru) ) (18) 
Equation (18) and the related integrals of (18) are what led us to believe that the Cole-Hopf transformation 
was quite general in nature. The particular form of (16) is 

U= — dr (19) 



where the diffusion equation is 

09 v d ( 86 



at rOr^O-r) (20) 
A simple example of the use of this method is to use the initial-value solution of the diffusion equation (a 
Dirac delta at t=r=0) 

0=-L e -*^ (21) 
2vt y ' 



Then the solution of (17) is found from (19) 



u=- t (22) 



This is the cylindrical version of Burgers' cartesian solution [4], u = x/t, which is the inviscid solution of 
(1), (solution 2.1, figure 5 in Benton [3]). Because (22) satisfies (17) in the limit that v — > (which is also 
the case of zero acceleration for fluid particles), this solution is the inviscid limit of the solution of Burgers' 
equation in cylindrical coordinates. Some solutions of (17), such as example B below, approach the sawtooth 
limit as t — > oo; but this is not surprising as it has been noted [2,13] that the sawtooth shape of dissipation 
layers is of fundamental importance in Burgers' representation of turbulence. 

It is extremely useful to derive solutions of the diffusion equation which are based on physical boundary 
conditions on the velocity field. To this end, equation (19) may be formally integrated to 

6(r, t) = k(t) exp J u{w, t) dw^j (23) 

where 

k(t)=0(O,t) (24) 
The velocity is to be initially specified over a given range so that the initial value of 9 is given by 

0(r,O) - o (r) = fco cxp J u (w) dw^j (25) 

for r < r < R. 

The general solution of the diffusion equation, (20), with azimuthal symmetry, is [19] 

/■OO />oo 

6(r,t)= / 9 {r')J {kr')Mkr)e- uk2t kr'dkdr' (26) 
Jo Jo 

The integral over k can be done once and for all as it is a special case of Weber's second exponential integral 
as discussed in Watson [21]. 
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cxp(-vtk 2 )J (kr)J (kr')kdk = -L exp ( J r l±J_ | j Q ( UL\ (27) 



2vt \ 4i/t ) \2vt J 
where the modified Bessel function of order zero may be written as 



/o(ar') = Jo(iar') 



(28) 
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which is a real function. Equation (26) may now be written as 



9{r,t) = 




4. Two Simple Examples 

We choose to work two relatively simple examples so as not to obscure the simple nature of the transform. It 
is surprisingly easy to pose apparently straight-forward problems that cannot be readily carried to completion 
because of extremely difficult quadratures. We will treat the more difficult case of the decay of an initially 
constant velocity wind in a future paper. 

The reader may have noticed by now that Burger's equation does not specify the density distribution except 
in a few limits besides the trivial one of constant density. 

A. Incompressible Flow 

The first case is that of incompressible flow where we take a source of water emanating from the origin of 
coordinates and spreading out in a circularly symmetric pattern. The initial flow is assumed to exist from 
very close to r=0 out to R, which is assumed to be larger than any other length scales in the problem. We 
then take the limit as R — > oo, in the same spirit as other studies of this initial- value problem on infinite 
intervals [2]. The speed must decrease as r~ 1 to conserve mass in an initially steady flow. In sum, 

M (r,0)=u o — (30) 



r < r < R 

so that (25) integrates to 

0o = fco(^) a (31) 

where 

a = — (32) 

Clearly, uq may be large at the small value ro in such a way that the product uo r o is well-defined. Equation 
(29) can now be integrated using Hankel's generalization of (27) as discussed in Watson [21], so that 



= W (l-f) (4*)-*E$* PS) 
where (|) is Pochhammer notation, 

(S). - S (S + 0(1 + s )~(f+- < 34 » 

and 



The full non- linear solution can now be derived from (19) and expressed in terms of the confluent hyperge- 
ometric function as 

af Af( fl /2 + l,2,-r>t) 
" (r '' )= 2t M(a/2,l,-rV4^) ^ 

or more clearly as 

/ 2 \ n— 1 

u r > *) = ™ 1 , 2\« — 366 

We see that the solution can be written as the product of the inviscid solution times an infinite sum of powers 
of the similarity variable ^ . We note the non-linear effect on higher powers of the similarity variable through 
the occurrence of the factor n in the numerator. The transient non-linearly steepens but diffusion dominates; 
note that the transient is not by itself a solution of Burgers' equation as we do not have the superposition 
theorem of linear theory. 

Figure 1 was created by non-dimensionalizing equation (17) in the usual fashion, i.e., r — > r / L,t — > i/t/l 2 , u — > 
u/U. Choosing U = uq and L = ro, then a is .5. This ratio of series is initially divergent and required 100 
terms in double precision on a Sun SPARCstation. The solutions are presented as isochrones, showing the 
evolution of the non-linear solution towards a steady-state. The boundary where the wind was turned off 
moves outward while diffusion smooths out the discontinuity. The wind initially drops off as 1/r so that 
the peak in the profile approaches the starting radius for the source as t — > 0. All curves approach the 
steady-state envelope of 1/r while diffusion flattens the leading edge. Asymptotically the velocity is zero 
everywhere. 

B. The Decay of An Initially Linearly Increasing Wind 

Example two is a wind emanating from the origin, again, but now the velocity linearly increases in magnitude 
out to some radius, r , chosen to be larger than any other length scale in the problem. We had in mind the 
initial theoretical treatments of radiatively-driven winds [6] from early-type stars, but our description is too 
simple to be applied to this problem. This is because of our assumption of cylindrical rather than spherical 
symmetry, and also because other effects are now known to be important in these winds [8]. 

The driving force is imagined to turn-off at t=0, and then we calculate the decay of the flow. In sum 

u{r,0)=u {) — = -P (37a) 
< r < r ; t = 

where 

to = — (376) 
u 

is half the time for the wind to flow from the source to the initial outer boundary. Both r , u are large, but 
t is not necessarily big. Then 

9 = ko e- br2 (38) 



Proceeding as before 

a k o 

= 7- cxp 

l + i 



(40) 



u(r,t) = — (41) 
We note that the solution approaches the Burgers' sawtooth shape as t — > oo. 

5. Conclusions 

We have given a relatively simple proof that Burgers' equation can be solved for curvilinear coordinate 
systems in three dimensions by assuming irrotational flow and using equation (16) with 9 derived from 
solutions of the diffusion equation. We solved two examples in cylindrical coordinates with circular symmetry 
to show how to implement the generalized solution. Our goal has been to extend the tools of workers who 
model the complicated interplay between non-linear steepening and viscous diffusion in waves. 
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Figure Caption 

The non- linear evolution of a wind with an initial spatial dependence of 1/r. The isochrones show the 
decay of the wind, and, especially, the diffusion at the leading edge where the wind was turned off as this 
discontinuity propagates outward. 



